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锅 积 累 对 艾 纳 香 内 生 菌 群落 结构 和 共 发 生 网 络 的 影 啊 


陈 娇 娇 ， 任 建国 ， 王 俊 丽 ” 
贵州 医科 大 学 公共 卫生 与 健康 学 院 ”环境 污染 与 疾病 监控 教育 部 重点 实验 室 ， 贵 阳 
561113) 

摘 要; 为 探究 器 官 锅 (Cd 积 累 对 艾 纳 香 内 生 菌 的 影响 ， 该 文采 用 168 rRNA RARA 
测序 技术 结合 分 子 生态 网 络 分 析 ， 研 究 不 同 外 源 锅 处 理 (0、2.0 mg'kg) 下 艾 纳 香 根 、 茎 、 
叶 Cd 积累 对 内 生 菌 群落 结构 特征 的 影响 。 结 果 表 明 : (1) 与 未 添加 外 源 Cd(0 mg:kg!, Cdo) 
相 比 ，Cd(2.0 mg*kg"*，Cd2) 处 理 促进 植株 生长 且 根 、 茎 、 叶 中 Cd 积累 量 为 叶 (16.75 mg kg 
1)» 25(11.99 mg'kgD> 根 (3.96 mg:kg-!))。(2) 内 生 菌 a 多 样 性 分 析 表 明 ， 各 器 官 内 生 菌 丰富 度 
(Sobs 指数 ，Ace 指数 和 Chao 指数 ) 和 多 样 性 (Shannon 指数 和 Simpson 指数 ) 均 以 根 最 高 ， 
茎 次 之 ， 叶 最 低 且 Cd2 处 理 各 嚣 官 内 生 菌 丰富 度 和 多 样 性 均 高 于 Cd0 处 理 。(3) 在 门 水 平 上 ， 
两 处 理 根 、 茎 、 叶 内 生 菌 均 以 变形 菌 门 、 放 线 菌 门 和 厚 壁 菌 门 为 优势 菌 群 ， 在 属 水 平 上 ， 
代 尔 夫 特 菌 是 Cd0 和 Cd2 处 理 各 器 官 的 主要 菌 属 ， 相 对 丰 度 分 别 为 53.0%~92.7% 和 
57.1%~89.2%; 艾 纳 香 根 、 茎 、 叶 内 生 菌 群落 结构 有 一 定 的 相似 性 ，Cd2 处 理 提高 根 、 苓 、 
叶 共 有 内 生 菌 属 的 比例 及 各 器 官 (根除 外 ) 独 有 内 生 菌 属 比 例 。(4)LDA Effect Size 分 析 表 明 ， 
相同 处 理 不 同 器 官 间 及 不 同 处 理 相 同 器 官 间 内 生 菌 属 存 在 差异 。(5) 宛 余 分 析 (RDA) 发 现 ， 
根 际 土壤 Cd. RES Cd 含量 与 内 生 菌 群落 结构 组 成 有 明显 的 相关 性 。(6) 共 发 生 网 络 分 析 结 
RER, Cd 积累 使 艾 纳 香 根 、 叶 内 生 菌 共 发 生 网 络 变 得 复杂 且 增 强 根 、 茎 物种 间 的 竞争 作 
用 和 叶 物 种 间 的 协同 共生 作用 。 综 上 所 述 ， 外 源 Cd 处 理 影响 艾 纳 香 各 器 官 内 生 菌 群落 结 
构 及 相互 作用 模式 。 
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Effects of cadmium accumulation on the structure and co- 
occurrence network of endophytic bacterial community in 


Blumea balsamifera 


CHEN Jiaojiao, REN Jianguo, WANG Junli* 

(School of Public Health, Key Laboratory of Environmental Pollution Monitoring and Disease 
Control, Ministry of Education, Guizhou Medical University, Guiyang 561113, China) 
Abstract: To probe into the impacts of organ cadmium accumulation on endophytic bacteria in 
Blumea balsamifera, the methods of high-throughput sequencing of 16S rRNA gene V3-V4 region 
and molecular ecological network analysis were employed to study the effects of Cd accumulation 
in roots, stems and leaves of B. balsamifera on the community characteristics of endophytic 
bacteria under different exogenous cadmium treatments (0 and 2.0 mg:kg'!). The results were as 
follows: (1) In comparison with the control group without exogenous cadmium addition (0 mg: kg: 
!, Cd0), the treatment group with cadmium spiking in soils (2.0 mg:kg'!, Cd2) promoted plant 
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growth and the cumulative Cd contents in root, stem and leaf, with the order of leaf (16.75 mg: kg; 
D > stem (11.99 mg:kg!) > root (3.96 mg: kg'!). (2) Alpha diversity analysis showed that the 
endophytic bacteria richness (Sobs, Ace and Chao indices) and diversity (Shannon and Simpson 
indices) for organs were the highest in roots, followed by stems and leaves under Cd0 and Cd2 
treatments. Additionally, the indices of richness and diversity of endophytic bacteria for each 
organ under Cd2 treatment were superior to those under Cd0. (3) At the phylum level, 
Proteobacteria, Actinobacteriota and Firmicutes were the dominant bacterial phyla in all the 
organs for both treatments; At the genus level, Delftia was the main bacterial genus with the 
relative abundance ranged from 53.096 to 92.7% and 57.196 to 89.296 in the plant organs of CdO 
and Cd2, respectively; Certain similarities existed among the endophytic bacterial community 
structures of roots, stems and leaves of B. balsamifera, and Cd2 increased the proportion of 
mutual endophytic bacterial genera in roots, stems and leaves and that of unique endophytic ones 
in each organ (except the root) (4) LDA Effect Size analysis showed that there existed the 
differences on endophytic bacterial genera residing in different organs within a group and also the 
same organ between groups. (5) Redundancy analysis (RDA) showed that the contents of 
rhizosphere soil Cd and plant organ Cd were significantly correlated with the composition of 
endophytic bacterial community. (6) Co-occurrence network analysis clarified that cadmium 
accumulation in B. balsamifera complicated the interaction network of endophytic bacteria 
occurring in root and leaf, and enhanced the competition among endophytic bacterial species in 
root and stem, and the symbiosis in leaf. In summary, exogenous Cd treatment affected the 
community structure and interaction mode of endophytic bacteria in B. balsamifera organs. 

Key words: Blumea balsamifera, endophytes, cadmium(Cd), 16S rRNA, co-occurrence 
network 


锅 (Cd) 是 我 国 农 田 重金 属 污染 元 素 之 一 ， 是 植物 非 必 需 元 素 。 随 着 工业 化 的 发 
展 ，Cd 通过 汽车 尾气 、 水 泥 厂 、 农 用 化 肥 和 冶炼 厂 粉尘 等 途径 排放 到 环境 中 ， 造 成 环境 
污染 (Zhou et al., 2017; Lin et al., 2022)。 Cd 污染 导致 植物 营养 缺乏 、 光 合 强度 降低 、 氧 化 
应 激发 生 等 ， 表 现 为 植物 生长 受 抑 制 ， 严 重 的 导致 植物 死亡 ( 张 星 雨 等 ，2021)。 中 药材 是 
中 医 临床 防 病 治 病 的 物质 基础 ， 在 临床 治疗 中 占有 重要 地 位 ， 其 质量 好 坏 直接 关系 着 临床 
疗效 。 已 有 研究 表明 ，Cd 在 多 种 中 药材 中 有 超标 的 风险 ， 且 对 中 草药 生长 发 育 及 活性 成 分 
积累 均 有 抑制 作用 ( 蒲 翔 等 ，2019;， 张 德 林 等 ，2019)。 

植物 内 生 菌 是 一 类 存在 于 植物 体内 、 不 对 寄主 植物 产生 伤害 的 微生物 类 群 (Hardoim et 
al., 2008)。 在 长 期 协同 进化 过 程 中 ， 内 生 菌 与 宿主 植物 形成 了 互利 共生 的 友好 关系 ， 表 现 
为 宿主 植物 为 内 生 菌 提供 稳定 的 生存 环境 和 营养 物质 ， 而 内 生 菌 则 通过 分 泌 代 谢 物 调节 宿 
主 植物 的 生理 活动 ， 进 而 促进 植物 从 土壤 环境 中 汲取 养分 ( 姚 雨 轩 ，2022)。 当 寄主 植物 遭 
遇 非 生物 (Cd 胁迫 时 ， 内 生 菌 主要 通过 激素 调节 (He et al., 2017)、 增 强 光合 作用 (Zhang et 
al., 2023)、 增 加 和 蛋白质 含量 (Zhou et al., 2021)、 增 强 抗 氧化 系统 (Kuramshina et al., 2018)、 强 
化 根系 细胞 壁 木质 素 沉 积 (Maslennikova et al., 2023) 以 及 与 宿主 协同 作用 (Xie et al., 2023) 等 
方式 缓解 重金 属 胁迫 ;或 通过 胞 内 积累 、 胞 外 固定 等 将 有 毒 重金 属 转化 成 低 毒 或 无 毒 形式 
(Long et al, 2021); 还 可 通过 分 泌 铁 载体 (Nakamoto et al., 2021)、 有 机 酸 (Li et al., 2017). 4& 
基 酸 (Latif et al., 2016) 和 提高 土壤 pH(Yuanli et al., 2022) 等 降低 重金 属 的 植物 有 效 性 ; 也 可 
通过 调节 植物 对 重金 属 吸收 、 转 运 基因 的 表达 ， 从 而 降低 或 增强 植物 对 重金 属 的 积累 
(Cheng et al., 2021; Qian et al., 2022). 

已 有 研究 表明 ， 超 积累 植物 不 同 器 官 或 组 织 内 细菌 群落 组 成 与 其 重金 属 积累 量 密切 相 
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1.6 高 通 量 数据 处 理 及 分 析 


值 20 以 下 的 碱 其 ， 人 允许 overlap 
OTU 聚 类 ， 去 除 所 有 样本 中 注释 到 的 叶绿体 和 线粒体 序列 ， 将 所 有 样本 序列 


纯化 PCR 产物 ， 测 序 工作 交 由 上 海 美 吉 和 4 


的 V3-V4 


和 806R :  5- 


区 进行 扩 增 ， 扩 增 体系 : 
) 0.8 nL， 下 游 引 物 (5 
因 组 DNA 10 ng， 补 


95 *C 30 min， 退 火 


使 用 fastp 软件 和 FLASH 软件 对 原始 涡 


区 最 大 错 本 


上 序 序列 进行 质 控 及 拼接 ， 过 滤 reads 尾部 质 
已 比率 为 0.2， 根 据 97% 的 相似 度 对 序列 进行 


利用 RDP classifier 分 类 网 站 和 Silva 16 rRNA 基因 数据 库 进 行 OTU 物 


EBE J 70%. iH% Alpha 多 


Ace 指数 )， 并 采用 Wilxocon 秩 和 检验 进行 Alpha 多 样 性 组 间 
(Linear discriminant analysis Effect Size) 确 定 同一 处 型 
f(LDA 值 >2)。 使 


属 水 平 差异 显著 的 内 生 


BEN: 
analysis, RDA) 探 讨 土壤 Cd 和 器 官 Cd 9] 3 483 8r 2 
genus 间 进 行 Spearman 相关 分 析 ， 得 到 相关 系数 和 矩阵 和 尸 值 矩 
Hochberg false discovery rate (FEDR) 方 法 矫正 上 述 计算 中 所 得 到 的 己 值 ， 根 据 genus 间 
Spearman 相关 系数 (>0.77 和 P<0.001 作为 阔 值 进行 物种 间 
软件 (0.10.1) 构 建 分 子 生态 网 络 (molecular ecological networks, MENS), 


1 AxyPrep DNA Gel Extraction Kit 
E 物 医药 科技 有 限 公司 完成 。 


三 | 
FH 


| 平 至 2 000, 


' 分 类 学 注释 ， 置 信 


样 性 (Shannon 指数 、Simpson 指数 、Sobs 指数 、Chao 指数 、 


连通 度 、 网 络 直 径 、 图 密度 、 


1.7 数据 处 理 


态 分 布 和 方差 齐 性 时 采 
2 结果 与 分 析 
2.4. 艾 纳 香 器 官 生物 量 及 Cd 


如 图 1: a 所 示 。 与 Cd0 处 到 


含量 


E 相 比 ，Cd2 处 理 下 


E 不 同 器 


差异 分 析 。 


] LEfSe 分 析 
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间或 不 同 处 理 组 间 相 同 器 官 


和 人 


m 


用 基于 


群 


ES 


平均 聚 类 系数 、 模 块 化 、 模 块 个 数 等 在 内 


使 用 SPSS 26.0 和 Excel 进行 数据 的 整理 分 析 和 绘图 


正 态 分 布 和 满足 方差 齐 性 时 ， 采 用 1 检验、 方差 分 析 (ANOVA) 法 
JIES% Wilcoxon 秩 和 检验 和 Kruskal-Wallis H SkRUVE US o 


欧 氏 距离 的 元 余 分 析 (redundancy 
落 结构 的 影响 。 


对 各 器 官 中 所 有 


并 采用 Benjamini and 


发 生 网 络 图 构建 。 利 用 Gephi 


并 计算 包括 网 


络 平均 


的 网 络 拓扑 属性 。 


-=> 


， 实 验 数据 


各 器 官 生物 量 均 增 加 且 


均 数 + 标准 差 。 服 从 
分 析 数 据 差异 ， 不 满足 正 


至 的 生物 


量 增加 明显 (P<0.05)，1 
纳 香 对 一 定量 的 Cd 有 了 耐 


受 性 。 


Cd 含量 测定 结果 (图 1: 


处 理 下 各 器 官 Cd 含量 显著 高 


于 Cd0 处 至 


此 可 见 2.0 mg:kg''Cd 处 理 对 艾 纳 香 


^E Keith 


MEHR, 


进而 说 明 艾 


b) 表 明 ， 艾 纳 香 各 器 官 Cd 积累 量 大 小 表现 为 叶 > 茎 > 根 且 Cd2 


E(P<0.05)。 此 外 ， 由 图 1: e 可 知 ， 艾 纳 香 对 外 源 


Cd 的 转运 系数 达到 15.82， 富 集 系 数 达 到 7.39， 说 明 艾 纳 香 对 Cd 具有 较 强 的 转运 和 富 集 
能 
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a. 器 官 生 物 量 ; b. 
水 平 。 


a. Biomass of the organs; b. Cd content of organs ; 


器 官 Cd 含量 ; 


c. 转移 系数 和 富 


Root Stem 


Leaf 


Cd0 and Cd2 represent 0 and 2.0mg: kg! Cd treatment levels, respectively. 


图 1 Cd 处 理 对 艾 纳 香 不 同 器 官 生 物 量 (a)、Cd 含量 (b)、 


Fig.l 


集 系数 ; Cd0 和 Cd2 分 别 代 表 0 和 2.0 mg'kg! Cd 处 理 


Translocation factor Bioaccumulation factor 


c. Translocation factor and bioaccumulation factor; 


富 集 系数 及 转运 系数 (c) 的 影响 


Effects of Cd treatment on biomass(a), Cd content(b), bioaccumulation factor and 


translocation factor(c) of different organs in Blumea balsamifera 


22 Cd 对 艾 纳 香 内 生 细 菌 Alpha 多 样 性 影响 


样本 


使 用 fastp. FLASH 软件 对 18 个 样本 的 原始 测序 结果 进行 
理 ， 共 获得 969 532 条 优化 序列 ， 其 中 Cdo 处 至 


过 滤 、 拼 接 、 算 选 等 后 续 处 
均 获 得 159 421 条 优化 序列 ， 


Cd2 


处 理 样本 平均 获得 149 711 条 优化 序列 。 在 97% 相 似 度 下 ， 序 列 被 注释 得 到 1 423 个 OTUs, 


Cdo 和 Cd2 处 理 分 别 得 到 616 个 和 807 个 OTUs( 表 
充分 地 反映 出 艾 纳 香 实验 样本 的 内 生 菌 


2 Br. WE 
艾 纳 香 根 、 


D. 2 il 


多 样 性 ， 本 实验 绘制 了 Chao 1 指数 的 稀释 


EX AM ERE 


本 的 测序 深度 能 否 


| 线 ， 如 


线 


内 生 菌 丰富 度 (Sobs 指数 、Ace 指数 、 


在 一 定 差异 ， 具 


体 表 现 为 根 中 内 生 菌 


直 增 高 后 呈 平 稳 状态 
Z. "AEN Alpha %7 M 


中 内 生 菌 多 样 性 


及 丰富 富 度 无 } 显著 差异 。 


为 根 > 葵 > 叶 。 
组 间 
高 于 Cd0 处 理 ， 


类 内 生 盏 以 使 其 适 


Alpha 多 样 怕 
特别 是 茎 中 内 生 菌 丰富 度 (P<0.05)， 
高 了 艾 纳 香 内 生 菌 Alpha 多 样 性 ， 也 进 


MERKE, X4 


应 体内 不 同 微 环 境 的 变化 。 


GR? 说 明 测 语 量 充足 ， 支持 后 续 的 分 析 。 
分 析 结 果 见 表 1。 相 同 处 理 条 件 下 ， 
Chao 指数 ) 和 多 样 


不 同 器 官 


性 (Shannon 指数 、Simpson 指数 ) 存 
丰富 度 及 多 样 性 明显 高 于 共和 叶 (CP<0.03)， 而 茎 、 叶 


内 香 各 器 官 内 生 菌 Alpha 多 样 性 呈现 


FE 结果 ( 表 1) 表 明 ，Cd2 处 理 下 艾 纳 香 各 器 官 内 生 菌 多 样 性 及 丰富 度 均 
由 此 可 见 外 源 Cd 添加 (2.0 mg-kg!) 提 


步 说 明了 宿主 器 官 Cd 积累 会 使 植株 招 


募 更 多 种 


X 1 Cd 积累 对 艾 纳 香 不 同 器 官 内 生 细菌 Alpha 多 样 性 影响 
Table 1 Effects of cadmium accumulation on Alpha diversity of endophytic bacteria in different 


organs of Blumea balsamifera 


序列 数 OTU 数 香农 指数 ”辛普森 指数 BK 
分 组 uH Sobs 指数 Ace 指数 Chao 指数 本 
Number of | Number 0n Nu» "M Shannon Simpson Coverage 
Group Organ Sobs indice Ace indice Chao indice NN ui 
Sequences of OTUs indice indice (%) 
Cd0 根 Root 47 481 297 434.67+70.67a 500.33+74.67a 487.71+75.39a 2.30+0.23a 0.31+0.69b 99.8 
Æ Stem 55 301 167 207.67+37.45b 232.21+45.09b 233.59+41.73b  0.58+0.11b C0.85+0.028a 99.9 
叶 Leaf 36 639 152 196.00+72.19b 216.76+80.84b 220.46+77.99b 0.53+0.12b 0.86+0.027a 99.9 
Cd2 根 Root 47 221 353 539.67+39.12a 599.35+67.51a 595.17+70.72a 2.58+0.15a 0.33+0.045b 99.8 
Æ Stem 52 800 265 428.33+83.68*ab 489.71+83.20*ab 490.13+72.67*a — 1.2740.43b 0.61+0.25ab 99.8 
叶 Leaf 49 690 187 276.67+106.98b 297.92+117.95b 279.92+112.70b 0.74+0.14b 0.80+0.054a 99.9 


T 


注 : * 代 表 不 同 处理 相 同 器 官 之 间 存 在 差异 ， 不 同 字母 代表 相同 处 理 不 同 器 官 之 间 存 在 差异 (P<0.05)。 


Noted: * represents the differences for the same organ between different treatments, and different letters represent 


the differences of different organs under the same treatment (P«0.05). 
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采用 渐 近 方法 评估 物种 丰富 度 。 
Species richness estimation by asymptotic approach. 


图 2 Chao 1 指数 的 稀释 曲线 图 


Fig.2 Rarefaction curve diagram of Chao 1 indice 


2.3 Cd 对 艾 纳 香 内 生 菌 群落 组 成 的 影响 

基于 门 水 平 的 内 生 菌 群落 相对 丰 度 分 析 表 明 ， 外 源 Cd 处 理会 影响 艾 纳 香 各 器 官 内 生 
菌 群落 组 成 (图 3: a)。 变 形 菌 门 (Proteobacteria)、 放 线 菌 门 (Actinobacteriota) 和 厚 壁 菌 门 
(Firmicutes) 等 在 各 器 官 中 占 绝对 数量 优势 。Cd2 处 理 根 、 茎 、 叶 中 放 线 菌 门 、 厚 壁 菌 门 、 
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175 个 (56.8%)。 
具有 一 定 的 相似 性 ， 但 不 同 处 到 
生 菌 对 各 器 官 不 同化 学 组 成 环境 适应 的 改变 。 此 外 ， 与 Cd0 处 理 相 比 ，Cd2 处 理 
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图 3 Cd0 和 Cd2 处 理 下 艾 纳 香 各 器 官 内 生 菌 群落 的 组 成 
Fig.3 Endophytic bacterial community in various organs of Blumea balsamifera under Cd0 and 
Cd2 treatments 
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图 4 Cdo 和 Cd2 处 理 下 艾 纳 香 各 器 官 内 生 菌 群落 构成 维 恩 图 ( 属 水 平 ) 


Fig.4 Venn diagrams of endophytic bacteria community in various organs of Blumea 
balsamifera under Cd0 and Cd2 treatments (genus level) 


通过 元 余 分 析 RDA) 识 别 内 生 菌 群落 与 根 际 土壤 Cd 和 器 官 Cd 含量 之 间 的 关系 。 由 
可知 ， 艾 纳 香 根 、 茎 、 叶 前 20 种 物种 中 ，Cd2 处 理 各 器 官 多 数 内 生 菌 丰 度 与 根 际 土壤 
Cd、 器 官 Cd 含量 呈正 相关 ， 包 括 根 中 芽孢 杆菌 属 、Mycobacterium、Sphingomonas、 
Pseudonocardia, Comamonadaceae, Tahibacter, RIR WREE, ZP Mycobacterium, 
Nocardioides, Actinoplanes, FWA WE. Curtobacterium, Bradyrhizobium, RMA Js ^5, 
叶 中 芽孢 杆菌 属 、Cowrtobactera ~ 8$ A EZ Pj J& ~ Acidovorax, Tahibacter 、 
Enterobacteriaceae, Massilia 等 ; 仅 有 少数 内 生 菌 与 之 呈 负 相关 ， 如 根 中 Allorhizobium- 
Neorhizobium-Pararhizobium-Rhizobium, EP RE~ Acidovorax, Massilia. WAIE W JE 
ZH Acidovorax, Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium, Massilia, KK K 
特 菌 属 ， 叶 中 Bradyrhizobium、 代 尔 夫 特 菌 属 、Rhodomicrobium、 Comamonadaceae « 
Nocardioides 等 ， 而 Cd0 处 理 各 器 官 内 生 菌 丰 度 与 根 际 土壤 Cd、 器 官 Cd 含量 的 关系 与 
Cd2 处 理 的 结果 相反 ， 上 述 结果 暗示 Cd 处 理 下 内 生 菌 组 成 和 结构 的 变化 对 艾 纳 香 Cd 富 集 
作用 有 着 重要 的 影响 。 
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(genus level) and Cd contents in Blumea balsamifera organs and rhizosphere soils 


24 Cd 对 艾 纳 香 内 生 菌 物种 差异 的 影响 
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一 处 理 组 内 不 同 器 官 间 及 不 同 处 理 组 相同 器 官 之 间 的 主要 差异 物种 。 结 果 发 现 ，Cd0 处 理 
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6 艾 纳 香 器 官 内 生 菌 响应 Cd 的 LEfSe 分 析 ( 组 内 ) 


Fig.6 LEfSe analysis on the response of endophytic bacteria to Cd in organs of Blumea 
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7 艾 纳 香 器 官 内 生 菌 响应 Cd 的 LEfSe 分 析 ( 组 间 ) 


Fig.7 LEfSe analysis on the response of endophytic bacteria to Cd in organs of Blumea 


balsamifera (Between groups) 


2.5 Cd 对 艾 纳 香 内 生 菌 共 发 生 网 络 的 影响 


共 发 生 网 络 总 体 拓扑 性 质 分 析 结 果 见 表 2。Cd0 AE PR. A. np 


114， 边 分 别 为 1 558、1 250、 
密度 分 别 为 0.119、0.184、0.209; 
140、152， 边 分 别 为 7200、998、 


别 为 162、117、 
23.632， 网 络 
别 为 254、 
41.303， 网 


发 生 网 络 规模 趋向 大 而 复杂 ， 茎 则 表现 


络 密度 分 别 为 0.112、0.103、0.274， 
出 相反 效果 。 


发 生 网 络 节点 分 
1347, TEENA 1.235: 21.368. 
Cd2 Ab FAR. zE. RRE NKT 4) 
3 139， 平 均 连 通 度 分 别 为 28.346、14.257、 
上 述 结果 表明 Cd2 AERE RENT 叶 内 生 菌 共 
Cd0 处 理 下 根 、 茎 、 叶 各 器 官 内 生 菌 


nea 


tk 发生 网 络 模块 性 数值 分 别 为 0.591、0.672、0.517， 模 块 数 分 别 为 29、15、14，Cd2 处 到 
F 各 器 官 共 发 生 网 络 模块 性 数值 分 别 为 0.654、0.693、0.346， 模 块 数 分 别 为 68、28、17， 
该 结果 表明 Cd2 处 理 下 根 、 茎 和 叶 共 发 生 网 络 内 生 菌 功能 单元 均 有 所 增加 ， 各 器 官 共 发 9 
网 络 模块 表现 (模块 性 ) 以 叶 网 络 模块 间 连 接 较 多 且 紧 密 ， 茎 次 之 。 此 外 ， 各 处 理 样 本 内 台 
菌 群 落 均 以 正 相 关 为 主 ， 但 随 着 Cd 的 施加 ， 根 、 茎 中 共 发 生 网 络 正 相关 连 线 比 例 减少 (人 
别 是 根 )， 而 叶 中 正 相关 连 线 比例 则 有 所 增加 ， 表 明 Cd 处 理 改变 了 各 器 官 内 生 菌 物种 间 的 
相互 作用 ， 具 体 表现 为 叶 中 协同 共生 关系 的 加 强 和 根 、 茎 中 物种 间 兖 争 关系 的 增强 。 

对 物种 进行 相关 性 网 络 互 作 分 析 ， 就 内 生 菌 共 发 生 网 络 结构 组 成 (图 8) 来 看 ，Cd0 处 理 
根 中 边 数 较 大 的 物种 为 Pseudoxanthobacter(42)、 Dokdonella(42)、 Novispirillum(42)、 
Bryobacter(42). Occallatibacter(42). Gemmatimonas(42). Azovibrio(42). Cryomorpha(42). 
Zymobacter(42) ~ Romboutsia(42) 等 ， 葵 中 边 数 较 大 的 物种 为 Planococcus(31) ~ 
Aquaspirillum(31). Turicibacter(31). Dactylosporangium(31). Devosia(31). Thauera(31). 
Haematobacter(31). Garicola(31). Agromyces(31). Altererythrobacter(31)5&, TW rp UBL 
的 物种 为 Planococcus(4l) .  Bdellovibrio(Al) . Dokdonella(41) ~  Ahniella(4l) ~ 
Hyphomicrobium(41) ~ Luteimonas(41) ~ | Cellulomonas(A1) ~ | Photobacterium(41) ~ 
Pseudarthrobacter(Al) ~ Reyranella(41) 等 ; Cd2 处 理 根 中 边 数 较 大 的 物种 为 
Sporacetigenium(81) ~ 乳 杆 W J& (81) ~  Novispirillum(81) ~ Truepera(81) ~ 
Altererythrobacter(81). Romboutsia(81). Reyranella(81). Gordonia(81). Brevundimonas(81). 
Dongia(81) 等 ， 茎 中 边 数 较 大 的 物种 为 Planifilum(28) 、  Luteimonas(28) ~ 
Dactylosporangium(28). Thauera(28). Streptococcus(28). Ohtaekwangia(28). Agromyces(28). 
Zymobacter(28) 、 Microvirga(28) ~ Rubrobacter(28) 等 ， 时 中 边 数 较 大 的 物种 为 
Deinococcus(70) . | Pseudoxanthobacter(70) ~ | Pedomicrobium(70) . | Planococcus(70) 、 
Rhodoplanes(70). Planifilum(70). Lechevalieria(70). Aquaspirillum(70). Nonomuraea(70) . 
Dokdonella(70) 等 。 由 此 可 见 器 官 Cd 积累 会 影响 各 器 官 内 生 菌 共 发 生 网 络 结构 和 物种 间 的 
互 作 关 系 。 
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表 2 艾 纳 香 各 器 官 共 发 生 网 络 特征 参数 


Table 2 Characteristic parameters of co-occurrence network of various organs of Blumea 


balsamifera 
Cd0 Cd2 
Network parameter 根 ES 叶 根 zE Hr 
Root Stem Leaf Root Stem Leaf 
T 162 117 114 254 140 152 
Node 
边 1558 1250 1347 7200 998 3139 
Link 
XE M MEABÉS HE 
平均 连通 度 19.235 21.368 23.632 28.346 14.257 41.303 
Average degree 
图 密度 
0.119 0.184 0.209 0.112 0.103 0.274 
Graph density 
模块 性 
- 0.591 0.672 0.517 0.654 0.693 0.346 
Modularity 
模块 数 
2 15 14 68 28 17 
Number ofmodules i 
正 相 关 1299 1097 1322 4212 838 2959 
Positive correlation (83.38%) (87.76%0) (90.7294) (58.595) (83.9794) (94.27%) 
负 相 关 259 153 125 2988 160 180 


Negative correlation (16.62%) (12.24%) (9.28%) (41.5%) (16.03%) (5.73%) 
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A. Cd0 处 至 
A. Cd0 treatment; 


ANE 
UM y 
E; B.Cd2 处 理 ; 红色 和 绿色 分 别 代表 正 相 关 和 负 相 关 (Spearman'”s 770.7. 和 Significant P«0.001). 


B. Cd2 treatment; 


Red and green links represent positive and negative correlation 


respectively(Spearman's r>0.7 and Significant P«0.001). 
图 8 Cd 积累 对 艾 纳 香 各 器 官 内 生 菌 群落 结构 的 影响 ( 属 水 平 ) 
Fig.8 Effects of Cd accumulation on community structure of endophytic bacteria (genus level) in 


organs of Blumea balsamifera 

3 Wit 
3.1 Cd 积累 对 艾 纳 香 器 官 内 生 菌 群落 结构 的 影响 

植物 内 生 菌 群落 易 受 非 生 物 和 生物 因子 所 调节 (Hallmann et al., 1997; Fuentes-Ramirez et 
al., 1999; Seghers et al., 2004; Liang et al., 2021)。 相 关 研 究 普 裔 认为 Cd 污染 导致 植物 内 生 
菌 多 样 性 下 降 ( 才 文 钰 ，2021; Li etal., 2023)， 但 在 黑 麦 草 (Lolium perenne) (Liang et al., 
2021)、 超 积累 型 东南 景 天 等 ( 邹 淑 华 等 ，2019) 相 关 研 究 中 发 现金 属 Cd 污染 导致 内 生 菌 多 
样 性 增加 ， 同 样 在 本 研究 中 ， 艾 纳 香 Cd 胁迫 积累 不 仅 能 促进 植株 生长 ， 且 提高 各 器 官 内 
生 细 菌 多 样 性 和 丰富 度 ， 该 现象 的 发 生 与 一 些 研 究 者 指出 的 内 生 菌 多 样 性 越 高 ， 宿 主 植物 
对 环境 胁迫 的 适宜 性 越 强 结果 相符 (Sinchez-L6pez et al., 2018)。 另 外 ， 本 研究 Cd0 和 Cd2 


处 理 艾 纳 香 根 中 内 生 细 菌 多 样 性 和 丰富 度 均 显著 高 于 茎 和 叶 ， 这 与 三 七 (Panax notoginseng) 
内 生 菌 器 官 分 布 结果 相似 (Liu et al., 2020)， 可 能 与 根 内 生 菌 主要 来 源 于 土壤 有 关 (Gao & Shi, 
2018)。 

本 研究 中 艾 纳 香 内 生 菌 群落 由 变形 菌 门 、 放 线 菌 门 和 厚 壁 菌 门 等 组 成 ， 其 中 变形 菌 门 
占 主要 地 位 ， 这 与 甘 松 (Nardostachys jatamansi) ( 李 莹 等 ，2022)、 冰 川 环 豆 (Oxytropzis 
glacialis) ( 许 国 琪 等 ，2021) 等 植物 内 生 菌 组 成 研究 结果 一 致 ， 说 明 植物 内 生 菌 在 较 大 分 类 
单元 下 基本 一 致 。 在 本 研究 中 ， 与 Cd0 处 理 相 比 ，Cd2 处 理 根 、 茎 、 叶 中 变形 菌 门 的 相对 
丰 度 降低 ， 而 放 线 菌 门 、 厚 壁 菌 门 、 酸 杆菌 门 的 相对 丰 度 升 高 ， 这 与 Liang 5:(2021). 425 
华 等 (2019) 等 研究 结果 一 致 。 而 从 属 水 平 上 来 看 ，Cd 积累 导致 艾 纳 香 根 、 茎 、 叶 中 优势 苗 
属 增 加 ， 如 根 中 的 代 尔 夫 特 菌 属 、 畏 氨 醇 菌 属 ， 茎 中 的 假 单 胞 菌 属 、 乳 杆菌 属 和 叶 中 的 芽 
TERT ER, RDA 分 析 结 果 也 表明 上 述 苗 属 与 器 官 Cd 含量 和 土壤 Cd 含量 的 密切 相关 性 。 
各 器 官 不 同 程度 的 Cd 积累 导致 内 生 菌 群落 结构 不 同 也 体现 在 东南 景 天 (Sedum alfredii) E(B 
淑 华 等 ，2019)， 如 超 积 累 型 东南 景 天 (HE) 葵 中 Alphaproteobacteria 、 叶 片 中 
Gammaproteobacteria、Negativicutes 和 Clostridia 4 个 纲 内 生 菌 的 相对 丰 度 随 Cd** 积 累 的 增 
多 显著 增加 。 与 根部 相 比 ， 超 积累 植物 Thlaspi caerulescens subsp. calaminaria 地 上 部 内 生 
茵 不 仅 有 了 耐 高 浓度 Cd*+、Zn** 等 特点 ， 还 存在 有 器 官 特定 优势 菌 属 ， 如 Methylobacterium 
(Lodewyckx et al., 2002)， 由 此 可 见 ， 内 生 菌 群落 结构 受 宿主 植物 和 土壤 环境 的 影响 。 

相关 研究 指出 ， 在 Cd 胁迫 下 ， 植 物 能 够 招募 特定 的 内 生 菌 种 群 ， 从 而 实现 高 Cd 抗 性 
和 植物 生长 促进 特性 ( 台 喜 生 等 ，2021; Liang et al., 2021)。 一 般 来 说 ， 招 募 的 特定 内 生 菌 
种 群 往往 具有 较 强 的 耐 Cd 性 ， 如 东南 景 天 从 根 际 土壤 中 招募 大 量 耐 Cd/Zn 内 生 菌 ， 包 括 
Burkpoideria、 芽 孢 杆菌 属 、Novospjpipgopiam 等 来 强化 其 对 重金 属 胁 迫 的 抵抗 (Wu et al., 
2020)。 本 文 研 究 结果 发 现 ，Cd 处 理 艾 纳 香 根 、 茎 、 叶 中 存在 大 量 的 差异 优势 物种 ， 如 芽 
孢 杆菌 属 、 假 单 胞 菌 属 、 乳 杆菌 属 等 。 相 关 研 究 指 出 ， 芽 孢 杆菌 属 能 够 强化 根系 细胞 壁 木 
质 素 的 沉积 ， 进 而 阻止 Cd 向 地 上 部 转移 ， 以 减轻 其 毒害 作用 (Maslennikova et al., 2023). 
假 单 胞 菌 属 可 通过 调节 东南 景 天 重金 属 ATPase(HMA)、 天 然 抗 性 相关 巨 哈 蛋白 (NRAMP)、 
ZRT/IRT-like 蛋白 (ZIP) 等 基因 的 表达 来 实现 对 Cd 转移 的 调控 和 耐性 (Chen et al., 2017)。 乳 
杆菌 属 的 胞 外 多 糖 通过 改变 水 稳 根 Cd 亚 细 胞 分 布 和 Cd 化 学 形态 形式 来 减少 Cd 向 地 上 部 
籽粒 的 转移 和 减轻 Cd 毒害 (Li et al., 2022)。 由 此 可 见 ， 艾 纳 香 各 器 官 中 更 多 差异 内 生 菌 可 
能 参与 艾 纳 香 生长 及 Cd2 的 吸收 、 积 累 和 解毒 生理 代谢 过 程 。 
32 Cd 积累 对 艾 纳 香 内 生 菌 共 发 生 网 络 结构 的 影响 
在 生态 系统 中 ， 微 生物 通常 以 复杂 的 网 络 结构 共存 和 相互 作用 ， 构 建 微生物 共 发 生 网 
络 为 探究 微生物 群落 内 潜在 的 相互 作用 提供 了 技术 支撑 (Chu et al., 2021)。 一 般 来 说 ， 网 络 
越 复杂 ， 和 群落 结构 越 稳定 (Li et al., 2022)。 艾 纳 香 各 器 官 内 生 菌 共 发 生 网 络 分 析 结 果 表 明 ， 
无 论 是 Cd0 或 Cd2 人 处理， 根部 内 生 菌 网 络 结构 都 比 茎 、 叶 复杂 且 物 种 间 以 正 相 关 为 主 ， 
与 Zhang 等 (2022) 研 究 结果 一 致 ， 即 在 两 种 农业 种 植 模式 中 ， 均 表现 为 番茄 根 中 内 生 旨 
网 络 比 茎 中 复杂 。 两 处 理 组 对 比 发 现 ，Cd2 处 理 下 根 、 叶 内 生 菌 共 发 生 网 络 比 Cdo 更 复杂 
(表现 为 更 多 的 边 数 和 较 大 的 平均 连通 度 和 图 密度 数值 )， 且 增加 根 、 茎 中 内 生 菌 共 发 生 网 
络 负 相 关 边 比例 和 叶 中 正 相 关 边 比例 ， 进 而 说 明 Cd 积累 导致 各 器 官 内 生 菌 之 间 的 相互 作 
用 关系 有 所 改变 ， 即 在 根 、 茎 中 物种 间 竞 争 关 系 的 增强 ， 而 在 叶 中 则 为 物种 间 共 生 关 系 的 
增强 ,但 从 总 体 而 言 ， 各 器 官 物种 间 的 共生 关系 是 占 主导 地 位 的 。 同 样 在 He 等 (2021) 研 究 
H, Ria Ate m S (Commelina communis) 内 生 菌 共 发 生 网 络 复杂 性 且 物 种 间 主 要 以 
正 相 关 关 系 存在 ， 而 超 积 累 植 物 刺槐 (Robinia pseudoacacia) 在 蚀 、 铜 、 铅 等 污染 状态 下 ， 
器 官 内 生 菌 共 发 生 网 络 趋 向 于 简单 ， 物 种 间 的 关系 并 不 以 正 相 关 作 用 为 主 ， 而 从 各 器 官 内 
生 菌 共 发 生 网 络 模块 特性 来 看 ， 叶 内 生 菌 共 发 生 网 络 比 根 更 复杂 ( 姚 雨 轩 ，2022)， 这 与 本 
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研究 中 Cd2 处 理 下 器 官 内 生 菌 共 发 生 网 络 特点 表现 基本 一 致 。 锅 污染 下 水 稻 器 官 内 生 菌 共 
发 生 网 络 分 析 表 明 ， 与 地 上 部 (模块 性 为 0.537) 相 比 ， 根 部 (模块 性 为 0.277) 积 累 更 多 Cd 的 
同时 其 内 生 菌 群落 结构 表现 为 更 复杂 化 (Zheng et al., 2022)， 这 与 Berry 和 Widder(2014) X 
述 的 复杂 网 络 能 更 好 应 对 环境 变化 观点 相 一 致 。 在 本 研究 中 ，Cd 积累 下 艾 纳 香 器 官 叶 的 模 
块 性 远 低 于 茎 和 根 ， 说 明 叶 部 内 生 菌 共 发 生 网 络 模块 间 联 系 紧密 且 不 同 物种 分 布 于 不 同 的 
模块 ， 有 利于 形成 特定 生态 位 相对 稳定 的 共生 内 生 菌 群落 以 适应 Cd 积累 胁迫 。 因 此 在 本 
研究 中 ，Cd 污染 可 能 通过 诱导 艾 纳 香 建 立 更 加 复杂 的 内 生 菌 共 发 生 网 络 和 增强 物种 间 的 正 
向 关联 强度 来 应 对 Cd 积累 的 影响 。 值 得 注意 的 是 ， 本 文 并 未 设置 Cd 对 艾 纳 香 生长 的 抑制 
浓度 ， 且 研究 样本 量 有 限 ， 存 在 不 足 之 人 处。 因此， 在 后 续 的 研究 中 将 进一步 改进 研究 方案 ， 
使 研究 变 得 更 加 深入 、 合 理 。 


4 结论 


(1) MiSeq 测序 表明 Cd 积累 提高 了 艾 纳 香 各 器 官 内 生 
多 样 性 总 体 表现 为 根 > 茎 > 叶 。 

(2) RDA 分 析 表 明 艾 纳 香 器 官 内 生 菌 群落 组 成 与 植株 体 器 官 Cd 含量 、 根 际 土壤 Cd 含量 
密切 相关 且 大 多 数 内 生 菌 与 植株 体 Cd 含量 、 根 际 土壤 Cd 含量 呈正 相关 。 

(3) 线性 判别 分 析 表 明 ， 艾 纳 香 根 中 内 生 菌 差异 物种 较 芭 和 叶 丰 富 。 此 外 ，Cd 积累 显著 
增加 根 、 茎 、 叶 中 Mycobacterium. Nocardioides. Bradyrhizobium, Sphingomonas. FUIF Té 
E~ Cellulomonas. 2£Í8EFFFAE. Mesorhizobium, Rhizobacter 等 属 的 耐 Cd 内 生 菌 数量 。 

(4) 共 发 生 网 络 分 析 表 明 ， 艾 纳 香 根 、 叶 中 Cd 积累 使 内 生 菌 共 发 生 网 络 变 得 更 加 复杂 。 
器 官 Cd 积累 增强 了 根 、 茎 物种 间 的 竞争 关系 和 叶 物 种 间 的 共生 关系 。 
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